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Abstract 


A combined finite element method/method of moments (FEM/ 

MoM) approach is used to analyze the electromagnetic scattering 
properties of a three-dimensional-cavity-backed apertu re in an infin ite 
ground plane. The FEM is used to formulate the fields inside the cav- 
ity, and the MoM ( with subdomain bases) in both spectral and spatial 
domains is used to formulate the fields above the ground plane. Fields 
in the aperture and the cavity are solved using a system of equations 
resulting from the combination of the FEM and the MoM. By virtue of 
the FEM, this combined approach is applicable to all arbitrarily 
shaped cavities with inhomogeneous material fillings, and because of 
the subdomain bases used in the MoM, the apertures can be of any 
arbitrary shape. This approach leads to a partly sparse and partly 
dense symmetric matrix, which is efficiently solved using a biconjugate 
gradient algorithm. Numerical results are presented to validate the 
analysis. 

1. Introduction 

The electromagnetic characterization of cavity-backed apertures is important in understanding their 
scattering properties and in electromagnetic penetration-coupling studies. Recently, there has been con- 
siderable interest in analyzing cavity-backed apertures in an infinite ground plane. Various analytical 
and numerical techniques have been applied for two-dimensional-cavity-backed apertures (refs. 1 to 6). 
For three-dimensional problems, mode matching has been used for rectangular- (ref. 7) and spherical- 
(ref. 8) cavity-backed apertures. Recently, a method of moments (MoM) modal (ref. 9) combined 
approach has been used to analyze apertures formed by rectangular cavities recessed in ground planes. 
These methods are restricted to cavities with regular shapes, where fields can be written in modal form. 
In reference 10, a boundary integral method is used to analyze the scattering from three-dimensional 
cavities by means of a connection scheme. Though useful for savings in computer memory, this method 
leads to dense matrices. Also, the accumulation of errors caused by the connection algorithm is not neg- 
ligible as the number of subsections increases. In the case of deep cavities, high-frequency techniques 
such as those proposed in references 11 and 12 could be effectively implemented. Unfortunately these 
techniques are not suitable when the cavity is filled with inhomogeneous materials. In reference 13 a 
finite difference time domain (FDTD) method is applied for large structures, but the FDTD method 
sometimes results in inaccurate results because of differencing, staircasing, and dispersion. Recently, 
Jin and Volakis (ref. 14) used a finite-element-boundary integral formulation that uses the boundary 
integral equation to formulate the fields external to the cavity accurately. 

In this paper, an approach similar to that described in reference 14 is followed. The volume of the 
cavity is discretized using tetrahedral elements, which automatically results in discretization of the aper- 
ture in triangular elements. The evaluation of the MoM integrals involves integration over triangular 
patches. These integrals are efficiently solved using the formulas given in reference 15. As an alterna- 
tive approach, these integrals are evaluated using the spectral-domain approach given in reference 16. 
Numerical studies indicate that both methods, when combined with the finite element method (FEM) 
give almost identical results, but the spectral-domain method becomes prohibitively time consuming for 
moderate to large aperture sizes. Comparisons of CPU times for both methods are given. Some exam- 
ples are considered with and without material fillings to validate the present analysis. Numerical results 
thus obtained are compared with the available data in the literature. 



Symbols 

A FEM/MoM matrix 

b excitation matrix 

ds surface integral 

dv volume integral 

E electric field 

E ; - incident electric field 

E z i z-component of E,- 

e unknown coefficient matrix 

ej unknown coefficients as defined in equation (18) 

F electric vector potential 

H magnetic field 

H, incident magnetic field 

H,. reflected magnetic field 

H scat scattered magnetic field 

H x j x-component of H, 

H y j v-componcnt of H ( 

H j z-component of H ; 

Hq j 0-component of H, 

H 0 j (^-component of H, 

H e s 0-component of H scat 

H ^ |, s (^-component of H scat 

j = -Pf 

k free-space wave number ~ 

A, 

k x Fourier transform variable with 

k y Fourier transform variable with 

k z complex propagation constant in z-direction 

Lj length of edge in tetrahedral element 

M magnetic current 

M x x-component of Fourier transform of M 

M y v-componcnt of Fourier transform of M 

h unit normal 

r position vector from the origin to the field point 

r' position vector from the origin to the source point 

r, 0, (|) spherical coordinates 

S surface 

S a aperture surface 

T vector testing function for volume elements 

T v vector testing function for surface elements 


respect to x 
respect to y 
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T x x-component of Fourier transform of T 

T y y-component of Fourier transform of T 

V volume 

W f vector basis functions for tetrahedral elements, i = 1 , 2, . . . , 6 

Wf vector basis functions for tetrahedral elements, j - 1 , 2, . . . , 6 

X, Y, Z rectangular coordinate axes 
x, y, z unit vector along X-, Y-, and Z-axis, respectively 

x bYb^b 
x 2’ Yb ^ 2 > 
x 3’ >3* Z3* 

x 4 , y 4 , z 4 coordinates of nodes of a tetrahedral element 
a polarization angle 

a m ,a n simplex coordinates associated with nodes m and n in tetrahedral element 
A integration over a triangular element 

A integration over a tetrahedral element 

relative permittivity 

r\ 0 intrinsic impedance of free space 

0,-, incident angles as shown in figure 1 

0 unit vector along 0-direction 

X free-space wavelength 

(X 0 permeability of free space 

(X,. relative permeability 

a scattering cross section 

a EE copolarized scattering cross section when E- polarized wave is transmitted 

a EH cross-polarized scattering cross section when E-polarizcd wave is transmitted 

c HE cross-polarized scattering cross section when H - polarized wave is transmitted 

<5 hh copolarized scattering cross section when H - polarized wave is transmitted 

(j) unit vector along ^-direction 

co angular frequency 

V del operator 


2. Analysis 

Figure 1 shows the geometry of the problem under consideration. The infinite ground plane is 
placed at z = 0- Consider that this geometry is illuminated by a harmonic plane wave (ref. 9) 


H, = 


(x/9 


, + yH y , + zH z i) e 


./' k , 


= + .) 


~J k i 


E / = THxk; 


(1) 

(2) 
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where 


k = ~k 0 (\ sin 0. cos <|). + y sin 0 ■ sin <|>. + z cos 0 ■) 


(3) 


H x j 


sin a cos 0 ■ cos A. + cos a sin A. 
— 




sin a cos 0. sin <h.- cos a cos A. 
— 


(4) 

(5) 


H :,i = 


-sin a sin 0. 


( 6 ) 


// 0 ■ = |H-| sin a (7) 

H^j = |H ; .| cos a (8) 

in which k a is the free-space wave number, r\ 0 is the tree-space intrinsic impedance, and a represents 
the polarization angle of the incident field. When a = 0, then H z ,• = 0, which corresponds to 
//-polarization, and when a = Jt/2, then E zi = 0, which corresponds to E-polarization. 

The electric field in the cavity caused by this incident wave satisfies the vector wave equation: 

Vxf— VxEVfc^E = 0 (9) 

VMy ) 


where p,. is the relative permeability and £,. is the relative permittivity of the medium. 

To solve the partial differential in equation (9) with the FEM, multiply equation (9) with a vector 
testing function T and integrate over the volume of the cavity. By applying suitable vector identities, we 
can write equation (9) in its weak form (ref. 17) as 


jjj -(VxT).(VxE)Jv-^e r ||jT.E dv 
v v 


JJJV.( T XJ.VX E )* 


( 10 ) 


Applying the divergence theorem to the integral on the right-hand side of equation (10) and using 
Maxwell’s equation V x E = -/Cop o p ( .H allows equation (10) to be rewritten as 

111 -( VxT )*( VxE )^ v - ^ £ r||| T,E dv = jwp 0 JJ(Txn).H ds (11) 

V V s 


where h is the outward unit normal to the surface enclosing the volume V and H is the total magnetic 
field at the surface. Because the tangential electric field is zero on the perfect-electric-conducting (PEC) 
walls of the cavity, the surface integral in equation (1 1) is nonzero only over the aperture opening in the 
infinite ground plane. In accordance with the equivalence principle, the cavity and the outside regions 
are decoupled by closing the aperture with a perfect conductor and introducing an equivalent magnetic 
current. By invoking image theory, H above the infinite ground plane consists of the incident field H ( -, 
the reflected field H,. from the ground plane, and the scattered magnetic field H scat from the equivalent 
magnetic current. Therefore, the total H field is 


H = H . + H ; . + H 


scat 


(12) 
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where, because of the perfectly conducting ground plane, H, = H r We can write equation (1 1) as 

111 rj~( V xT )*( v xE )^v - « E dv - y (0|0. o J' J (T X n) • H scat ds = 2;cop o J J (T x h) • H,. ds (13) 

V r V s a s a 

where S indicates the integration over the aperture surface. 

At this point, the problem can be divided into two parts, the first being the FEM part involving the 
discretization and evaluation of volume integrals on the left-hand side of equation (13) and the second 
being the MoM paid involving the discretization and evaluation of the surface integrals. We can write 
equation (13) as 


7 v + I s = 


(14) 


where 


I v = HI -J-(VxT).(VxE)dv-^e r |J|T.Edv (15) 

v v 

r s = -/®F 0 1 1 (X x n)*H scat ds (16) 

I e = 2ywp o ||(Txn).H / ^ (17) 


2.1. FEM Formulation 

To discretize the volume integrals in equation (13), the volume of the cavity is subdivided into 
small-volume tetrahedral elements. The electric field is expressed in terms of the edge vector basis func- 
tions, which enforce the divergenceless condition of the electric field implicitly (ref. 17): 

6 

E = £ e l W i (18) 

( = l 

where e ■ are the unknown coefficients associated with each edge of the tetrahedral element and 

W = LX a Voc -a Va ,) (19) 

i z v m n n v 7 


where m and n are the nodes connected to form edge i, L ■ is the length of the edge, and a ]u and a n are 
the simplex coordinates associated with nodes m and n. The simplex coordinates for the nodes of a tetra- 
hedron element are given by Silvester and Ferrari (ref. 1 8) as 


V 

V 

V 
V 4 

“4 = y 


a, = 


a 0 = 


a, = 


( 20 ) 

( 21 ) 

( 22 ) 

( 23 ) 
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where V is the volume of the tetrahedron given by 

1 x i y i z i 
1 x 2 y 2 z 2 

1 x 3 y 3 z 3 

1 x 4 v 4 z 4 

and V] , V 2 , V 3 , and V 4 are given by 

1 x y z 
1 x 2 y 2 z 2 
1 x 3 y 3 z 3 

1 x 4 v 4 z 4 

1 x l >T Z 1 

V = I 1 x -V Z 

6 1 x 3 y 3 z 3 

i ^ 4 y 4 z 4 

1 x i - v t z i 
1 X 2 y 2 Z 2 

1 x y z 
i * 4 y 4 z 4 

1 X 1 - V 1 M 
1 X 2 y 2 Z 2 

1 x 3 y 3 z 3 
1 x y z 

The testing function T is taken to be the same as an element of equation (18) that is, 





( 24 ) 


(25) 


(26) 


(27) 


(28) 


T = W j (j =1,2 6) (29) 

Substituting equations (18) and (29) into equation (15) yields 

6 

7 v = rr L lll (Vx W/)*(VxW 7 .) e . < /v-^e r JJJ(W i .W J .)« f rfv O' =1,2 6) (30) 

^ ; = l a A 

where ▲ indicates the integration over the volume of a tetrahedral element. 

The evaluation of these integrals over a tetrahedral element is given in reference 17. The resulting 
element submatrices are assembled over all the elements in the cavity volume to form a symmetric, 
sparse matrix of the FEM. 
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2.2. MoM Formulation 

The discretization of a cavity volume into tetrahedral elements automatically results in the discreti- 
zation of the aperture into triangular elements. In accordance with the equivalence principle (ref. 19), 
the fields in the two regions can be decoupled by closing the aperture with a perfect electric conductor 
and introducing the equivalent magnetic current 


M = E x z (at z = 0) (31) 

over the extent of the aperture. Hence, from equation (18), the magnetic current can be expressed in 
terms of the unknown coefficients associated with the tetrahedral elements as 

3 

M = 2> ( .(W.xz)| z = 0 (32) 

i = 1 


Similarly, 


T v = T x z = W j x z (j = 1, 2, 3) (33) 

The scattered magnetic field due to the magnetic current given in equation (31) above the infinite 
ground plane is given by 


H 


scat 


— — (Vx VxF) 


(34) 


where 


F = 7T dx ' dy ' <35) 

where r is the field point and r' is the source point with coordinates x and y'. In writing equations (34) 
and (35), the principle of image theory is used. 

2.3. Spectral-Domain Approach 

Writing equation (35) in terms of plane waves (ref. 19) and substituting equation (32) into 
equation (35) follows a procedure similar to that in reference 16 and allows equation (16) to be written 
as 


l = 


_ J_ 


471 ' 


r r (K — k x ) - ~ r r 

JJ - JL jf L M x (-k x , -k y )T x (k r k y )dk x dk y - J J 


(k k ) - 

-k y )T x (k x , k y )dk x dk y 
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(k 2 -k 2 J. - f f (k k ) - 

y -M,(~k x , -k y )T y (k x , k y )dk x dk y -J J — ir^~M x (-k r , -K)T y (k x , k v )dk y dk 


k ^ y 

K z 


Z 


(36) 


where 


k z = 


kl-kl-k 2 y for k] > (k 2 x + ky) 


-jjk 2 x + ky-k 2 0 for k 2 0 <{k 2 x + k 2 y ) 


(37) 
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and (M x ,M y ) and (T x , T y ) are the x- and y-componcnts of the Fourier transforms of M and T. The 
Fourier transform of a function F(x,y) is defined as 


rr j(k r x + k v y) 

F (k x , k y) = j \F{x,y)e x y ' dxdy (38) 

A 

Equation (36) is evaluated over each triangular element, and contributions from all the elements 
over the aperture are assembled to form a full matrix. This matrix is added to the FEM matrix to form 
the left-hand-side coefficient matrix. The numerical implementation of equation (36) is done by obtain- 
ing closed-form expressions for the Fourier transforms of the basis function over a triangular element 
(ref. 20). 

2.4. Spatial-Domain Approach 

The scattered magnetic field in equation (34) can be rewritten as (ref. 19) 

H scat = -T7-[^ F + V ( V * F )] < 39) 

Substituting equation (39) into equation (16) gives 

I s = -k 2 0 JJt v .F^-|Jt s .V(V.F) ds (40) 

The integrals in the above equation are evaluated over all the triangles on the aperture surface using 
the procedure outlined in Rao et al. (ref. 21). When calculating the mutual term over the triangles, a 
13-point Gaussian quadrature formula developed for triangles (ref. 22) is used. While evaluating the self 
term of these integrals, closed-form expressions given by Wilton et al. (ref. 15) are used. The dense 
matrix thus formed is assembled over all the triangles and combined with the FEM matrix using a global 
numbering system. 

2.5. Excitation Vector 

With the testing function defined in equation (33), the excitation integral in equation (16) can be 
evaluated by substituting for H, and integrating over each triangular element of the aperture. The 
element contributions are assembled over all the elements of the aperture to form the excitation vector. 

2.6. Solution Algorithm 

By adding the contributions from the FEM and MoM integrals and evaluating the excitation vector, 
equation (14) can be written in matrix form as 

Ae = b (41) 

where A is a symmetric matrix that is partly sparse (because of the FEM) and partly dense (because of 
the MoM). The column vector e is the unknown coefficient vector to be solved. Matrix b is the excita- 
tion vector. Because of the sparsity and symmetry of matrix A, only the nonzero elements of half of the 
matrix (upper half or lower half, including the diagonal) are stored to minimize the memory demand. A 
biconjugate gradient algorithm (ref. 23) solves equation (40) efficiently. 
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2.7. Scattered Cross Section 

Once the electric field E is found and hence the magnetic current M on the aperture, the far zone 
scattered field can be computed. 


H «C)I 


it 


K o e 


'l f , 2ttr 


J’J'(00 + $$).M(x',/) 


, , jk„ sin9(x' cos <Sf + y' sin 0) , 


dx' dy' (42) 


where (r, 0, (f)) are the usual spherical coordinates of the observation point. The scattering cross section 
is then given by 


2 H scat (r) 

a = lim 47tr~ (43) 

r -»“ |H f (r)| 

where 

|H„(r)| 2 = |H e ,,| 2 +|« t ,J 2 (44) 

|H,(r)| 2 = |»e. i | 2 + |ff*, j | 2 (45) 

The total scattering cross section is defined by equation (43); however, in most measurements either 
E-polarizcd or //-polarized waves are transmitted, and the E- polarized and //-polarized scattered far 
fields are measured separately. In such cases, the scattering cross section may be defined as 


c hh = li m 47tr 


Z\ H ty,s 


r — > o° o 

H <\ ),i 


a he = lim 471 r 


2| H e, s 


r — > oo TT 


r 4 2 \H e , s 

= lim 47tr 

I H 0 ; | 


a eh - lim 47tr 


! e, i 


2 1 


r ^“ H Q,i 


(46) 


(47) 


(48) 


(49) 


3. Numerical Results 

Backscatter cross section calculations were performed for various cavities to validate the analysis 
presented in this paper. These numerical results were compared with those available in the literature 
(refs. 9, 10, and 14). A commercial software package, COSMOS/M (ref. 24), models and generates the 
mesh for the geometries used for numerical computations in this paper. A discretization size of 0. 1A, 
generally suffices to achieve convergent results. Figure 2 shows the geometry of the rectangular cavity 
used in the numerical calculations shown in figures 3 to 6. 

To compare the spatial- and spectral-domain methods used to evaluate the MoM integrals, the back- 
scatter cross section of an empty cavity with dimensions 0.7a. x 0. 1 A x 1.73A for incident angle 0/ = 40° 
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was calculated as a function of (f) using both methods. The results are presented in figure 3. Both meth- 
ods lead to almost identical results, but the spectral-domain method took more time to fill the MoM 
matrix than the spatial-domain method did. In this example, a total of 1234 tetrahedral elements were 
used to discretize the cavity, which resulted in 20 triangular patches at the cavity aperture. The spectral- 
domain approach required 840 sec of CPU time to fill the matrix even after using the symmetry property 
of the MoM matrix, whereas the spatial-domain approach required only 5 sec of CPU time to fill the 
MoM matrix after using the symmetry property of the MoM matrix. The above calculations were done 
on a CONVEX C220 computer. Because of the large difference in CPU times required to fill the MoM 
matrix, the spatial-domain approach was used to carry out the rest of the numerical calculations. The 
spectral-domain approach became prohibitively time consuming as the number of triangular elements in 
the aperture increased. Backscatter patterns of this cavity for H- and E-polarizations were also com- 
pared with those results given in reference 14 and are shown in figure 4. 

Backscatter cross section was also calculated for a long and narrow cavity (a = 2.5X, b = 0.25 A,, and 
c = 0.25 A) and is given in figure 5. There were 944 tetrahedral elements used to discretize the cavity 
volume. This discretization procedure resulted in 154 triangular elements in the aperture. Comparison 
with the results given in reference 9 shows very good agreement. 

Figure 6 presents the backscatter cross section for a cavity with dimensions a = 1.0A, b = 0.25A, and 
c = 0.25A, and filled with material having £,. = 7.0 — j 1.5 and p r = 1.8 —jO.l. For this cavity, 2861 
tetrahedral elements were used in its discretization. This discretization procedure resulted in 200 trian- 
gular elements in the aperture area. In this study, a discretization size of 0. 1 k/j£ r generally resulted in 
convergent results for dielectric -filled cavities. The data computed with the present method agree well 
with those reported in reference 14. Backscatter cross-section calculations were performed for a rectan- 
gular cavity (fig. 7(a)) partly filled with a lossy dielectric material. The numerical data calculated were 
compared with those of reference 10 and are presented in figure 7(b). 

The following examples prove the flexibility of the present method. Figure 8(a) shows the geometry 
of a circular cavity with a circular aperture. Backscatter cross-section calculations for this cavity were 
made with the present method and are presented in figure 8(b). The data obtained with the present 
method agree well with those given in reference 10. As a second example, an air- filled rectangular cav- 
ity with a circular aperture was considered. The geometry of the problem is shown in figure 9(a). Back- 
scatter cross-section calculations for E-polarized and //-polarized incident fields were performed and 
are presented in figure 9(b). The numerical results obtained with the present method are compared with 
unpublished results obtained with a combined mode matching/MoM method. Because of the discontinu- 
ity of the circular aperture a discretization size of 0.03 5 A. was used in this example to achieve conver- 
gent results. A third example illustrates the flexibility of the present method for cavities filled with 
multidielectric materials. A rectangular cavity with a circular cylinder of dielectric material in its center 
and surrounded by a different material was considered. The geometry of the problem is shown in 
figure 10(a), and a discretization size of 0.05A, was used. Backscatter cross-section calculations were 
performed for both E- and //-polarizations of incident fields and are presented in figure 10(b). 

4. Conclusions 

A combined finite element method/method of moments (FEM/MoM) technique can be used to ana- 
lyze the electromagnetic scattering properties of a three-dimensional-cavity-backed aperture in an infi- 
nite ground plane. Both spectral- and spatial-domain approaches can be used to calculate the MoM 
integrals, but the spectral-domain approach is prohibitively time consuming for large apertures. The 
spatial-domain approach is much faster and as accurate as the spectral-domain approach. This approach 
is capable of calculating backscatter cross sections for cavities with arbitrary shapes and material fill- 
ings and with arbitrarily shaped aperture openings. The symmetry and sparsity of the matrix to be 
solved are completely exploited for savings in computer memory, and the biconjugate gradient 
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algorithm is used to save computer time. Numerical results obtained with this method show good agree- 
ment with the available data in the literature, proving its validity. 
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Figure 2. Geometry of rectangular-cavity-backed aperture. 


13 




Spatial-domain method, ^ j.j: 




Figure 3. Backscatter cross section of air-filled cavity versus incidence angle for a = 0.7X, b = 0.1 K, 
c = 1.73X, and 0, = 40°. 
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FEM/MoM results, o g j; 
FEM/MoM results, & uj 



(a) ^-polarization. 



Figure 4. Backscatter cross section of air-filled cavity versus incidence angle for a = 0.7X, b = 0.1 K, 
c = 1.73X, and 0, = 40°. 
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(a) Cavity geometry. 



Figure 8. Backscatter cross section of air-filled circular - cavity. 
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(b) Backscatter for E- and //-polarization. 

Figure 9. Backscatter cross section of air-filled rectangular cavity with circular aperture a = 0.5X, 
b = 0.5X, c = 0.3A, and radius 0.2A,. 


20 




(a) Cavity geometry. 



Figure 10. Backscatter cross section of inhomogeneously filled cavity with a = 0.5k, b = 0.5k, c = 0.3k, 
£,.| = 2.0, e,. 2 = 1.5, and inner dielectric radius 0.2k. 
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